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ABSTRACT 

We investigate the role of Alfven waves in the core-collapse supernova (SN) explosion. We assume 
that Alfven waves are generated by convections inside a proto-neutron star (PNS) and emitted from 
its surface. Then these waves propagate outwards, dissipate via nonlinear processes, and heat up 
matter around a stalled prompt shock. To quantitatively assess the importance of this process for 
the revival of the stalled shock, we perform ID time-dependent hydrodynamical simulations, taking 
into account the heating via the dissipation of Alfven waves that propagate radially outwards along 
open flux tubes. We show that the shock revival occurs if the surface field strength is larger than 
~ 2 x 10 15 G and if the amplitude of velocity fluctuation at the PNS surface is larger than ~ 20% 
of the local sound speed. Interestingly, the Alfven wave mechanism is self-regulating in the sense 
that the explosion energy is not very sensitive to the surface field strength and initial amplitude of 
Alfven waves as long as they are larger than the threshold values given above. 
Subject headings: MHD - supernovae: general - waves 
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1. INTRODUCTION 

The most promising scenario of collapse-driven super- 
nova (SN) is currently supposed to be the so-called de- 
layed explosion by neutrino heating (jKotake et alJl2006t 
I Janka et alj 120071 and references therein): the prompt 
shock, which was generated by core bounce and stalled by 
neutrino emissions and dissociations of nuclei, is heated 
up and revived by neutrinos coming out of the proto- 
neutron star (PNS), leading eventually to a SN explosion. 
Under spherical symmetry, however, no successful explo- 
sion has been obtained so far even though up-to-date 
micro-physics, such as equation of state an d weak inter- 
actio n rates, have been fully incorporated (iJanka et al.l 
2007). Various other effects have also been explored 
over the years. The implications of stellar rotation 
and different sorts of hydr odynamical instabili ties have 
bee n extensively studied. jKotake et al. (2006 |) and see 
also IMarek fc Jankal ()2007t ): [Mezzacappa et al.l (|2007f ) for 
very recent progresses.) 

Magnetic fields are drawing much attention 
of researchers these days although the his- 
tory of research is q uite long. After some 

pioneerin g papers (iLeBlanc fc Wilson! 
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had been forgotten for a while because it was realized 
that very strong magnetic fields are required to affect 
the supernova dynamics, which was supposed to be 
unrealistic at that time. The situation changed with the 
observational evide nce of strongly magnetized neutron 
stars or magnetars ((Thompson fc Duncan |[T996f l and the 
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progress in theoretical under standing of the magne to- 
rotational instability or MRI (jBalbus fc Hawlevl[l99ll ). 

The possible importance of MRI in the gener- 
ation of magnetic fi eld in core-co ll apse SNe was 
first pointed out by I Akivama et al.l (l2003ft. A lot 

of papers have been published (lYamada fc Sawail 

2004 [Thompson. Quataert. fc Burrows I " 120051 : 
Moiseenko. Bisnovatvi-Kogan. fc Ardeljanl 

120061 : iBurrows et ail 120071 : 
12007( 1 to investigate possible 
roles of magnetic fields in the supernova dynamics, 
especially on their kinematical aspects such as magnetic 
pressure and torque induced by rapid rotations i n super- 
nova cores. iThompson. Quataert. fc Burrows I {2005), 
on the other hand, considered a MHD turbulence that is 
possibly induced by MRI as a source of viscosity to tap 
free energies stored in differential rotations. Although 
the rapid rotation of stellar core is prerequisite in these 
studies, it may not be so easy to obtain in the presence 
of magnetic fields according to recent stellar evolution 
models because the transfer of angular momentum is 
efficient and the core may rot ate rather slowly just prior 
to th e gravitational collapse ((Heger. Woosley. fc Spruitl 
l2005h 5 . 

More recently, yet another supernova mechanism was 
put forward: sound waves generated by PNS oscil- 
lations of mainly g-mode nature, which are prob- 
ably induced by turbulence caused by the stand- 
ing accretion shock instability (SASI), heat up mat- 
ter through nonlinear dissipations, revive the stalled 
shock wave and prod uce explosions at very late times 
(|Burrows et al.l [2006H . Since this acoustic mechanism 
consists of several steps, whose efficiencies appear to 
be not very high (Yoshida. Ohnish iTfc Yamadal 120071 : 



IMarek fc~J anka 2007), the viability of the mechanism is 
still controversial and further explorations from various 



5 It is worth noting that Blondin & Mezzacapoa (2007) recently 
pointed out a new mechanism to generate the rotation of PNSs. 
They claimed that the non-axisymmetric SASI might be a source 
of the angular momentum of PNSs. If true, no rotation may be 
needed prior to the collapse to account for the spin of neutron stars. 
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view points are needed. 

At present the supernova mechanism is still elusive in 
spite of these extensive efforts. The above-mentioned 
research trend naturally leads us to the exploration of 
still another type of waves that may also contribute to 
the supernova explosion: Alfven waves. If the super- 
nova core is magnetized, the oscillations of PNS will 
emit not only sound waves but also Alfven waves. It 
is also possible that Alfven waves are excited by con- 
vections in PNS, w hich have been demonstrated to exist 
after core bounce (jKeil. Janka. fc Miiller I [l996h . This 
is quite analogous to what h appens in the Sun (e.g. 
iSuzuki fc Inutsuka \\2(M l2(M hereafter SI05; SI06). A 
fraction of these waves propagate outwards and dissi- 
pate later through nonlinear processes, which then will 
heat up matter and may lead to the revival of the 
stalled shock. Such a scenario as a supernova mech- 
anism has not been studied on a quantitative basis so 
far, altho ugh it has been discu s sed in a qualitative man- 
ner (e.g. I Wheeler et all 120001 : IWooslev fc Jankal [20051 : 
iBurrows et al.ll2006f ) or i n the context of the nucleo syn- 
thesis in ^-driven winds ([Suzuki fc Nagataki1l2005h . It 
should be emphasized that it is the dissipations of fluc- 
tuating components of magnetic fields associated with 
Alfven waves that heat up matter and revive the shock 
wave, whereas in most of the previous papers, torques 
and pressures exerted by total magnetic fields were the 
key players. Since the main aim in this paper is not 
to construct a realistic model but to elucidate the char- 
acteristics of the Alfven wave heating in a quantitative 
manner, we employ ID simplified but dynamical simula- 
tions. 

2. MODELS 

In this section, we describe the scenario we have 
in mind and explain nonlinear dissipation processes of 
Alfven waves. Based on these we give the basic equa- 
tions and employed approximations in detail. 

2.1. Basic Picture 

Convections around the neutrinosphere that roughly 
coincides with the surface of a PNS are induced by neg- 
ative gradients of lepton-fraction and entropy about a 
few tens of milliseconds after core bounce. According 
to hydrodynamical simulatio ns of PNS convections by 
iKeil. Janka. fc Miiller I (j!996f ). the amplitude of velocity 
fluctuations near the PNS surface is ~ 4 x 10 3 km s _1 
on average and becomes as high as 10 4 km s -1 . Re- 
cent simulations also reported that SASI in accreting 
flows mig ht excite a comparable order of su rface fluc- 
tuations (jOhnishi. Kotake. fc Yamada I [2006) . And not 
to mention, the non-radial oscillations of PNS are an ob- 
vious possibility as a source of Alfven waves. 

These activities in PNS generate various modes of 
waves with the amplitude of velocity fluctuations men- 
tioned above, which then emanate from its surface and 
propagate outwards. If the magnetic field is sufficiently 
strong, two types of magnetic waves are important: 
Alfven wave and fast wave. The Alfven waves prop- 
agate along field lines and are less subject to dissipa- 
tions thanks to their incompressive character. On the 
other hand, the fast waves propagate almost isotrop- 
ically and can traverse field lines although they suffer 
more damping owing to the compressive nature. There- 



fore, Alfven waves can transport energy farther away 
along magnetic filed lines unless they form closed loops. 
The fast waves could become important, however, in the 
equatorial region if PNS rotates rapidly and the field 
lines are tightly wound up. In this paper, we assume 
that open-field regions prevail in the post-bounce core 
SN core and study the energy deposition by the dissipa- 
tions of Alfven waves through nonlinear processes and its 
implications for the shock revival. 

The energy source in this mechanism is eventually the 
gravitational energy that is released by the matter ac- 
cretion onto PNS. Some of the released gravitational en- 
ergy is converted to the kinetic energy of surface fluc- 
tuations in the PNS by convections, SASI and other in- 
stabilities. Then the surface fluctuations excite various 
waves in PNS, some of which are Alfven waves traveling 
outwards along field lines. This sequence of processes is 
quite similar to what Burrows et al. (2006) proposed in 
their acoustic mechanism. 

As will be shown later, we assume strong magnetic 
fields of the magnetar scale, ~ 10 15 G, on the PNS sur- 
face. We do not specify the origin of the strong mag- 
netic fields in this paper (but see ^4.11 for related dis- 
cussions.) and treat the field strength as a free param- 
eter. As mentioned above, we focus on the propagation 
of Alfven waves and pay attention only to magnetic flux 
tubes that are extended beyond the stalled shock be- 
cause the Alfven wave carry the energy parallel to a field 
line; if the magnetic field is closed near the PNS surface, 
Alfven waves cannot transfer energy to the stalled shock. 

Since, as a first step, the aim of the paper is not to 
construct a realistic model, which should be the next 
step, but to understand the essential features of the 
Alfven wave heating in the post bounce core, we make 
the models as simple as possible, incorporating only min- 
imum ingredients. Since stellar rotation is not an indis- 
pensable player in the mechanism considered here, we 
just neglect it and assume radial magnetic fields as the 
simplest configuration for open field lines in this paper. 
The field strength is given by 

r 2 

B r = Brfi—^, (1) 

where B r ,o is the field strength at the PNS surface, 
r = r a . We perform one-dimensional (ID) simulations, 
ignoring the effects of neighboring magnetic fields. 

Even if the stellar rotation is fast, our models will be 
still applicable locally to the polar region, where the cen- 
trifugal force is minimum, although the simple prescrip- 
tion for the expansion factor of flux tubes, oc r 2 , adopted 
in this paper may need further elaboration. In the equa- 
torial region, on the other hand, the validity of the as- 
sumption in this paper depends not only on the rotation 
period but also on the reconnection efficiency among the 
field lines frozen into the accreted matter. This issue will 
be discussed again in £14.11 

2.2. Dissipation of Alfven Waves 

The effects of the Alfven waves on accreting matter are 
twofold: (1) the heating of matter by the Alfven wave 
dissipations and (2) the extra pressure exerted by the 
Alfven waves. We take into account only the former in 
this paper. The Alfven wave with linear amplitude is 
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non-dissipative due to the incompressive nature. How- 
ever, nonlinear Alfven waves suffer various dissipation 
processes. The nonlinearity, w, of Alfven waves can be 
defined as 

_ Sv± 5B ± 

where <5vj_ and 5B± are the amplitudes of velocity and 
magnetic field, the subscript '_L' denotes the transverse 
directions with respect to B r , and va = B r /^/4ivp is 
the Alfven speed. Here we have used the relation, 
6v ± = SBx/^/Airp, that the A lfven wave satisfies (§10 
of lLamers fc Cassinelli Ill999| ). Even though the initial 
amplitude at the launch from the PNS surface is in the 
linear regime, w -C 1, it grows thanks to the decrease of 
B r (Equation JlJ). Eventually, the Alfven wave becomes 
nonlinear, w ~ 1, and various dissipation processes set 
in. 

The excitation of compressive waves by nonlinear 
mode conversions is a route of the dissipation. If 
the Alfven wave is not strictly circularly polarized, 
the magnetic pressure, 5B^/8ir, fluctuates along with 
B r and induces longitudinal compressive motions, most 
of which correspond to slow magnetohydrodynamical 
(MHD) waves (SI05;SI06). Even if the Alfven wave is cir- 
cularly polarized, it is subject to the parametric decay in- 
stability, which generat es outgoing slow MHD waves and 
incom ing Alfven waves (Goldste irilll978tlTerasawa et al.l 
1986). The velocity amplitudes of the slow waves are 
also amplified as they propagate outwards and the den- 
sity decreases. Eventually the wave fronts steepen to 
form shocks and heat up matter. 

If B r has a transverse gradient (along _L direc- 
tion), fast MHD waves that propagate perpendic- 
ularly to B r are also excited from Alfven waves 
(|Nakariakov. Roberts, fc Murawski I Il997f) . These fast 
waves also heat ambient matter by the shock formation. 
Moreover, the tr ansverse inhomogeneity o f B r leads to 
the phase mixing (|Hevvaerts fc P riest 1983), which is an- 
other channel of the dissipation of Alfven waves. 

The turbulent cascade may also work in the dissipa- 
tion of Alfven waves. In the PNS, the Alfven speed 
is not constant along the radial magnetic field because 
of the changes in both B r and density. Then, incom- 
ing Alfven waves are excited from the outgoing ones 
by the deformation of wave shapes and the interac- 
tions between the outgoing and incoming Alfven waves 
lead to the formation of smaller scale (large wave num- 
ber) structures mainly in t he perpendicular directions 
(jGoldreich fc Sridhar 1 Il995). This turbulent cascade 
to higher wave number proceeds up to the dissipation 
range, where resistivities (for magnetic field fluctuations) 
and/or viscosities (for velocity fluctuations) become im- 
portant. Then, the energy that Alfven waves are carrying 
is finally transferred to the ambient matter. 

As a result of these various dissipation processes, the 
nonlinearity, w, of the outgoing Alfven waves is saturated 
at a certain level. According to dynamical simulations 
of solar (SI05;SI06) and stellar ([Suzuki 112001 winds, the 
saturation level is found to be w < 0.3—1, which we apply 
in this paper to the Alfven waves in the supernova core. 
It should be noted that the saturation with a constant 
w implies that SB± itself decreases as the Alfven wave 
propagate further outwards, since B r is declining (see 



Equation {5])). This means the dissipation of the energy 
of Alfven waves and results in the heating of the matter. 

2.3. Formulation 

We evaluate the Alfven wave heating by considering 
the conservation of wave energy under the WKB, or short 
wave length, approximation. This treatment enables us 
to obtain the estimations in a simple manner without 
solving the fully nonlinear MHD equatio ns. The conser- 
vation of wave energy can be expressed ([Jacques 111977ft 
as 

BE 

— f + V • F w - v ■ VP W = -pq w , (3) 
at 

where p, v are the density and velocity of the accretion 
flow and £ w is the wave energy density, F w is the wave 
energy flux, P w is the wave pressure, and q w is the dis- 
sipation rate of wave energy per unit mass. As for the 
Alfven wave, £ w = p{8v 2 x ) = (SBj_)/4-7r and P w = £ w /2, 
where the bracket (■ ■ •) stands for the average over a pe- 
riod of the Alfven wave. The last term on the left hand 
side denotes the rate of work done by the Alfven waves 
on the accretion flow, which we neglect in this paper as 
stated above. The right hand side represents the energy 
deposition by the Alfven waves to the accreting matter. 
When q w > 0, the matter is heated up by the wave dis- 
sipations, whereas the Alfven waves travel without dissi- 
pation if q w = 0. 

It is convenient for later use to introduce an adiabatic 
const ant, the so-call ed wave action, H w , defined as fol- 
lows ([Jacques II1977D : 



dP w 
dr 



va 



VA 



(4) 



We neglect relativistic corrections because they are minor 
outside the PNS. It should be noted that the wave action, 
-Hw, instead of _F W , is conserve d in moving m edia. The 
specific form of H w is given as ([Jacques Ill977[ ) 



(SBl) (v A + v)(v A + v) 



4tt 



''a 



(5) 



In this paper we assume the steady propagations of 
Alfven waves and neglect the time derivative, Jj, in 
Equation ([3]). This is valid when the Alfven transit time 
is shorter than a typical time-scale of the system, which 
will be discussed in $4.2\ 

The initial amplitude of Alfven waves at the PNS sur- 
face is supposed to be an order of convective velocities at 
the surface of PNS, which are suggested to be a fraction 
of the sound speed 6 : 

Svx,o = ec S:0 , (6) 

where c s .o is the sound speed at the PNS surface and 
we choose e — 0.1 — 0.3 in our simulations. Because 
c s , is «10% of the light speed, e = 0.1 — 0.3 corre- 
sponds to 5v± t o — (0-3 — 1) x 10 4 km s _1 , which is 
comparable to the values obtained in hy drodynamical 
simulations ([Keil. Janka. fc Miiller I [l996h . For a suffi- 
ciently large background magnetic field at the PNS sur- 
face, B t q > 5 x 10 14 G, the initial wave amplitude is 

6 This is very similar to the surface convection in the Sun. The 
observed granulation speed is 1 — 2 km s — 1 , while the sound speed 
at the photosphere is ss 5 km s _1 . 
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small in the sense that w < 0.1. In such a condition, 
the Alfven wave travels outwards without dissipation 
(q w = 0) near the PNS surface and its amplitude evolves 
according to the following relation: 
(6Bl)(v A + v r ) 2 



TABLE 1 
Summary of Simulations. 



,or Q , 



(7) 



47T Ua 

where i? w ,o is the wave action at its surface, and v r is 
the radial velocity of the background accretion flow. 

As the Alfven waves travel outwards, the nonlinear- 
ity, in, increases because of the expansion of the radial 
magnetic flux tube (Equation {T])). The dissipation of 
Alfven waves eventually sets in when they reach the non- 
linear regime as discussed above and the nonlinearity of 
the Alfven waves is saturated at a certain level: 

w w a. (8) 
In this paper we adopt a constant a = 0.5 as a standard 
saturation level, based on our previous results on the 
solar and stellar winds (SI05;SI06;Suzuki 2007), which 
showed that w(< 1) is more or less constan t or very 
slowl y varying as a function of r (see also e.g. iHollwegl 
119731 for steady-state modeling). This simple prescrip- 
tion of the constant saturation level is expected to in- 
corporate phenomenologically all the complex physical 
processes of Alfven wave dissipations discussed in §2.21 
and provide us with a reasonable heating rate. 

For the steady state, the energy dissipation rate, q w , 
can be obtained from the conservation of the wave action 
without referring to the details of the nonlinear processes 
that are responsible for the dissipations. Using the re- 
lation in the dissipation region, SB± oc B r oc (ro/r) 2 
(Equation {!])), the heating rate, q w (erg g~ 1 s _1 ), can be 
expressed as 

va a 2 B r d 2 n — 



va + v r inp dr 
An important point here is that we can evaluate q w from 
the local distributions of p, v r , and B r . 

2.4. Simulations 

Incorporating the above formula for the heating by 
Alfven waves, we perform ID time-dependent simula- 
tions of the post-bounce evolution s of SN core. We 
use th e 15M Q progenitor star of IWooslev fc Weaver! 
(fl995l) as an initial conditio n, and employ a numerical 
code (Sumiyoshi et al. 2005) to solve general relativistic 
hydrodynamics and neutrino transport, adding the extra 
heating term given by Equation ([9]) in the energy equa- 
tion. The Alfven wave heating (<j w ) is switched on at 
100 ms after core bounce. We have in mind that this 
is the time for the develop ment of (magneto-)convection 
that drives Alfven w a ves (|K cil, Janka, & Muller 119961 : 
lAkivama et ail [200l iMasada. Sano. fc Takabe I 120061 ) . 
Incidentally, we define the PNS surface as the position of 
the density, po — 10 11 g cm -3 , which, as a consequence, 
becomes smaller as the PNS contracts. 

We have three parameters, -B r ,0j an d ot in the above 
prescription. We fix the saturation level, a — 0.5, which 
controls the location of wave dissipations. The other two, 
the field strength, B r ^, at the PNS surface and the nor- 
malized initial amplitude, e, of perturbations, determine 
the wave energy injected from the PNS surface. In the 
following, we explore the condition on these two param- 
eters for the shock revival by simulating nine models in 
Table [Q 



Model 


B r>0 (G) 


e 


Explosion 


^cxp 




M cu t 


I 


1 x 10 15 


0.1 


No 








II 


1 x 10 15 


0.2 


No 








111 


1 x 10 15 


0.3 


No 








IV 


2 x 10 15 


0.1 


Marginal 








V 


2 x 10 15 


0.2 


Yes 


1.2 


0.08 


1.38 


VI 


2 x 10 15 


0.3 


Yes 


1.6 


0.10 


1.37 


VII 


3 x 10 15 


0.1 


Yes 


0.33 


0.04 


1.41 


VIII 


3 x 10 15 


0.2 


Yes 


1.5 


0.10 


1.38 


IX 


3 x 10 15 


0.3 


Yes 


2.2 


0.14 


1.38 


Note. — The explosion energy, ejecta mass and PNS mass 
are denoted by E exp , M e j and M cu t, respectively. The unit 
of -Eexp is 10 51 erg and the units of M c j and JWcut are Mq. 



3. RESULTS 

Figure Q] displays the success or failure of shock re- 
vival in the plane of the field strength at the PNS sur- 
face, -B r ,0i and the initial amplitude of velocity pertur- 
bation, e. Table [1] gives more detailed information on 
the ejected mass, M e j, explosion energy, E exp , and rem- 
nant (« PNS) mass, M cut , for the models that obtain the 
shock revival under the current approximation. They are 
estimated at the time t — 200 ms after core bounce as 
follows. We first define the ejecta as the collection of 
mass shells with positive total energy, which is the sum 
of kinetic, internal and gravitational energies. Then M e j 
and E exp are obtained as the sums of mass and total en- 
ergy, respectively, of each mass shell that comprises the 
ejecta. In so doing, the non-relativistic expression is em- 
ployed for the energy. It is clear that an explosion with 
-Ecxp > 10 51 erg is obtained if the magnetic field at the 
PNS surface is strong, B r _ > 2 x 10 15 G, and if the initial 
amplitude of Alfven waves is larger than a certain value, 
e(= Sv±/c s ) > 0.2. Note that the surface field strength 
is of the same order as those inferred for magnetars. 

Figure presents a typical Alfven wave-driven shock 
revival, in which the result for model V (red lines) 
is superimposed on the original result for the non- 
magnetized spherically symmetric model (black lines) 
(|Sumivoshi et al.ll2005l ). The trajectories of mass shells 
are plotted against the time from core bounce. The 
figure clearly demonstrates the shock revival by the 
Alfven wave he ating for the otherwise failed neutrino 
heating model (|Sumivoshi et al.1 I2005D like those dis- 
cussed in many previous papers (e.g. iKotake et alJl2006L 
and references therein). 

Figure [3] shows the rates of Alfven wave heating and 
neutrino heating (top panel) along with the velocity 
distribution (bottom panel) at t = 100 ms. The top 
panel demonstrates that the Alfven wave heating oper- 
ates mainly in the vicinity of the stalled shock wave and 
dominates the neutrino heating. The main reason for the 
localization of the Alfven wave heating in the Eulearian 
frame is the trapping of Alfven waves. The propagation 
speed of the outgoing Alfven wave is VA+v r in this frame. 
It rapidly decreases from 2 x 10 4 km s _1 at r = 100 km to 
ps km s _1 at r = 300 km for model V (bottom panel of 
Figure[3|), for example. The Alfven waves cannot travel 
further outwards and are trapped inside r < 300 km in 
model V, so that they spend a long time there to damp al- 
most completely. In the Lagrangian frame that moves at 
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Fig. 1. — Status of each model. The circles and crosses re- 
spectively corresponds to the models that result in shock revival 
and the models that do not give shock revival. The triangle is a 
marginal case. The radii of the circles in the successful cases are 
scaled by -E cxp . 

the inflow velocity of the accreting matter, on the other 
hand, the Alfven waves propagate outwards at the speed 
va- This means that the inflowing matter is heated up 
by the dissipation of Alfven waves just when it reaches 
the vicinity of the stalled shock wave. 

The luminosity, La, of Alfven waves at the PNS sur- 
face can be estimated under the assumption of spherical 
symmetry as follows: 

La = ^Po^Lo^A.o^r 2 , 




This implies that the emission of Alfven waves for ~ 100 
ms gives the energy injection of ~ 10 51 erg. In most cases, 
almost all the energy of Alfven waves is absorbed thanks 
to the trapping of Alfven waves just mentioned. This 
is also confirmed by the comparison with the simula- 
tion results. We can calculate the total heating rate (erg 
g _1 ) by multiplying the heating rate per unit mass (erg 
g -1 s _1 ; the top panel of Figure[3]) by the density (~ 10 9 g 
cm ) and volume of the heating region, 47rr 2 Ar, where 
Ar(« 100km in Figure [3]) is the thickness of the heat- 
ing region. We thus obtain the total heating rate of 
ps 10 52 erg s _1 for model V, which indicates that the in- 
jected Alfven wave luminosity is mostly used for heating 
up the stalled shock. 

10* e J 
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FlG. 2. — Time evolutions of mass shells. Model V with 
Alfven wave heating (B r ,o = 2 X 10 15 G Sz e = 0.2; red lines) is 
superimposed on the original non- magnetized model (black lines). 
The time is measured from the core bounce. The dashed lines show 
the locations of shock front. 
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FlG. 3. — Top: The heating (cooling for negative values) rates 
from neutrino (dot-dashed) and Alfven waves (solid) of model V 
at t = 100ms. The arrow indicates the location of the shock front. 
Bottom: The distributions of — w r (> 0) (solid), v\ (dashed), and 
D r + v\ (dot-dashed) at t = 100ms for model V. 

The regions with fast accretion velocities are preferen- 
tially heated up and eventually start to move outwards, 
provided La is sufficiently large. Once the stagnated 
shock wave is re-launched, the heating is reduced because 
Alfven waves become untrapped again. For larger La the 
shock revival occurs earlier and the duration of heating, 
Ata, is shorter. As a result, E cxp (roughly oc LaAta) is 
not very sensitive to B r o and e. In fact, although La of 
model IX is larger than that of model V by more than 
a factor of 3, the difference in E exp is less than a factor 
of 2. In this sense the Alfven wave mechanism is self- 
regulating. Interestingly, the acou stic wave mechanism 
is also claimed to be self-regulating ([Burrows et al.ll2006f ) 
though the regulating mechanism is different: the gener- 
ation of the acoustic waves continues until the shock is 
revived and matter ceases to accrete. 

Model VII is exceptional among the explosion cases, 
giving a very small explosion energy. The Alfven wave 
heating operates in a much outer region in this case be- 
cause Alfven waves become nonlinear, SB±/ B r > a, only 
after crossing the shock wave owing to the large B ri0 and 
small e. A sizable fraction of La leaks out of the stalled 
shock wave and a tiny amount is ejected with a quite 
small -Ecxp- 

The models with B r ^ = 1 x 10 15 G produce no ex- 
plosion. This is first because La itself is small (mainly 
models I & II) owing to the small B r .o and second be- 
cause the dissipation of Alfven waves occurs too early 
(mainly model III). Note in particular that La of model 
III is larger than those of the explosion cases V and VII. 
In this case Alfven waves dissipate in the inner region 
and the temperature is increased there. As a result, the 
energy deposited by Alfven waves is mostly converted to 
neutrino emission in this case. 

4. DISCUSSION 

In this section, we discuss more in detail the valid- 
ity of the assumptions and approximations employed for 
the background magnetic field and the formulation of 
Alfven wave propagations in this paper. 

4.1. Magnetic Field 

Since we have seen that strong magnetic fields (£? r .o > 
2 x 10 15 G), which are of the magnetar scale and much 
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larger than those for the ordinary radio pulsars, are 
necessary to revive the stalled prompt shock by the 
Alfven wave heating, one may think that the mechanism 
considered in this paper is only applicable to magnetar- 
forming supcrnovae and the Alfven wave is not a major 
ingredient in the ordinary SN explosion. We cautiously 
note, however, that this strong magnetic field is required 
not for the ordinary neutron stars as we observe them 
but for the PNSs in their very infancy. As a matter of 
fact, there is a speculation that the magnetic fields of 
nascent PNS may be temporarily very strong and then 
decrease to the 'normal' value (~ 10 12 G) by energy re- 
leases occurring during the SN exp losion and later evo- 
lution (| Wheeler. Meier, fe Wilson Il2002t) . If this is true, 
the Alfven wave mechanism may work in a larger popu- 
lation of core-collapse SNe. 

The origin of such strong magnetic fields is still contro- 
versial. One possibility is referred to as the fossil origin 
hypothesis: the strong magnetic field in compact stars is 
simply a consequence of the compression of the magnetic 
field that already exists in OB progenitors prior to the 
gravitational collapse. In fact, several magnetic massive 
stars have been observed to have a magneti c field whose 
avera ge dipole-field stren gth is ~ 1000G dNeiner et al.1 
l2003t iHubrig et all 120061 ; iDonati et all I2006D . The to- 
tal magnetic flux of these stars is comparable to that o f 
a typical magnetar (jFerrario fe Wickramasinghe I l2006h . 
which implies that an additional generation and/or am- 
plification of magnetic fields will not be necessary to 
obtain a highly magnetized compact remnant for these 
stars. 

Another possibility is an ampl ification of weak 
magnetic fields by the MRI (jAkivama et al.l 120031 : 
iMasada. Sano. fc Takabel [2006) , in which stellar rota- 
tion plays a key role, winding poloidal fields and driving 
the instability. In this scenario, toroidal magnetic fields 
are efficiently produced, whereas in the case of the fossil 
origin, we expect the radial component of magnetic field 
is dominant. 

Since the Alfven wave carries energy along a field line, 
we are interested only in open magnetic flux tubes that 
are extend beyond the stalled shock wave. As the sim- 
plest configuration, we have considered radial magnetic 
fields and neglected the toroidal component in this pa- 
per. As mentioned above, the approximation is justified 
if the magnetic field is of fossil origin and the progenitor 
is a slow rotator. Even if the progenitor core is a rapid 
rotator, our models will be still applicable to the polar 
region, where the effect of rotation is not strong and the 
toroidal magnetic fields are less important, although the 
radial dependence of the field strength may need elabo- 
ration. 

In the equatorial region of a rapidly rotating super- 
nova core, on the other hand, the situation is much more 
complicated. Field lines are not directed radially in gen- 
eral and some of them may be closed, as expected for 
the dipole configuration. In addition, the continuous 
downward advection of magnetic fields may cause recon- 
nections in the PNS, which in turn will open up some 
field lines again. In any case the toroidal component will 
be do minant over the radial component ( {Burrows et al.l 
l2007t ) and we need to include the effects of these spiral 
magnetic fields as well as rotation itself in discussing the 
Alfven wave heating in supernova cores quantitatively, 



which will be the future task. 

4.2. Alfven Wave Propagation 

The treatment of Alfven waves in this paper is ad- 
mittedly a crude approximation. We employ the non- 
relativistic, steady-state, and WKB approximations for 
describing the propagation of Alfven waves. Among 
the assumptions, the non-relativity is adequate for the 
Alfven waves launched from the PNS surface since the 
relativistic corrections are indeed minor outside the PNS. 
The steady-state approximation is also reasonable be- 
cause the Alfven transit time is shorter than the ex- 
pansion time of accreting matter; while the expansion 
time-scale of the ejecta is ~ 50 — 100 ms, the time for 
the Alfven wave to travel from the PNS surface (r w 50 
km) to the wave trapping region around the stalled shock 
(r « 200 km) is - 10 ms for d a «2x 10 4 km s" 1 . 

The WKB approximation is acceptable if the wave- 
length is shorter than the scale height of the background. 
The typical period of Alfven waves generated in the 
PNS is supposed to be t ~ 1 ms, corresponding to the 
dynamical time-scale. Then, the wavelength becomes, 
A « v a t ~ 2 x 10 4 km s" 1 x 10~ 3 s - 20 km. Since the 
density scale height is shortest near the PNS surface and 
is Hp w 30km, Alfven waves might be partially reflected 
there due to the deformation of the wave shape. In more 
detailed study, this effect should be taken into account 
together with the wave pressure ignored in this paper. 

5. CONCLUSION 

In this paper we have studied the matter heating by 
Alfven waves in the post-bounce supernova core and its 
implications for the shock revival. In order to elucidate 
the essential features of the mechanism quantitatively, we 
have done a couple of ID dynamical simulations, neglect- 
ing rotation and toroidal magnetic fields but employing 
the Alfven-wave heating rate based on our model of the 
nonlinear damping of Alfven waves. 

We have found that if the surface magnetic field 
strength is > 2 x 10 15 G and if the surface velocity fluctu- 
ation is e > 0.2, which corresponds to Sv±fi >6x 10 3 km 
s _1 , the stalled shock is revived by the Alfven wave heat- 
ing with a canonical explosion energy, E cxp > 10 51 erg. 
The current mechanism is self-regulating in the sense 
that the explosion energy is not very sensitive to the sur- 
face field strength and initial velocity fluctuation as long 
as they satisfy the above conditions. The above strong 
magnetic field is not a requirement for the ordinary neu- 
tron stars as observed but for the PNSs in their infancy. 
If magnetic fields decay through their subsequent evo- 
lution (Wheeler et al.2002), the Alfven waves may play 
an important role in a larger population of core-collapse 
SNe. 

It has been also found that the wave trapping is essen- 
tial in localizing the Alfven wave heating in the vicinity 
of the stalled shock wave as well as in regulating the ex- 
plosion energy. In fact, if the magnetic field is weaker 
(< 2 x 10 15 G), the Alfven wave heating takes place 
much closer to the PNS because the Alfven wave becomes 
nonlinear earlier on at smaller radii. Then no shock re- 
vival occurs because the dissipated energy is mostly lost 
by neutrino cooling; neutrino emissions are enhanced in 
this case. If the initial velocity fluctuation is smaller 
(e < 0.2) with a stronger magnetic field, >3x 10 15 G, on 
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the other hand, most of Alfven waves propagate through 
the stalled shock region with only a small amount of en- 
ergy being deposited, which then results in a weak ex- 
plosion with E cxp < 10 51 erg. It is noted, however, that 
even in these cases the Alfven wave heating will be still 
important in supplementing the neutrino heating. 

As a first step, the models presented in this paper 
have much room for sophistication. Among other things, 
as mentioned repeatedly, rotation and toroidal magnetic 
fields should be somehow taken into account. In real- 
ity, the Alfven wave mechanism probably works in co- 
operation with these processes. This will require multi- 
dimensional numerical modellings and will be the future 
work. 
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